Non-Equilibrium Kondo Model with Voltage Bias in a Magnetic Field 



P. Fritsch^'Q and S. Kehreini 

' Physics Department, Arnold Sommerfeld Center for Theoretical Physics, and Center for NanoScience, 
Ludwig-Maximilians-Umversitdt, Theresienstrasse 37, 80333 Munich, Germany 

(Dated: March 16, 2009) 

We derive a consistent 2-loop scaling picture for a Kondo dot in both equilibrium and non- 
equilibrium situations using the flow equation method. Our analysis incorporates the important 
decoherence effects from both thermal and non-equilibrium noise in a common setting. We calculate 
the spin-spin correlation function, the T-matrix, and the magnetization as functions of applied 
magnetic field, dc-voltage bias and temperature. In all these quantities we observe characteristic 
non-equilibrium features for a nonvanishing external voltage bias like Kondo splitting and strongly 
enhanced logarithmic corrections. 



I. INTRODUCTION 
A. Motivation 

The Kondo effect was first observed in the 1930's while 
measuring the resistivity of "pure" metals. Upon lower- 
ing the temperature one finds a minimum in the resis- 
tivity of nonmagnetic metals containing a small concen- 
tration of magnetic impurities. When lowering the tem- 
perature even further the resistivity increases and satu- 
rates at a finite value at zero temperature. Systematic 
experimental and theoretical analysis showed that this 
effect is due to a screening of the impurity spin by res- 
onant scattering of conduction band electrons leading to 
an enhanced electron density around the impurities. By- 
passing electrons scatter off these so called spin compen- 
sation clouds leading to an enhancement of the resistiv- 
ity. The Kondo model has become a paradigm model 
for stron g-co upling impurity physics in condensed mat- 
ter theory^. It has been solved exactly using the Bethe 
Ansatz'^ ^, however, dynamical quantities like the impu- 
rity spectral function are not easily accessible within this 
framework. Many other numerical and analytical meth- 
ods have been developed sin ce that can get around this 

limitatioiPanmiiiniiiiiiEni. 

Experiments on quantum dots in the Coulomb 
blocka de regim e have revived the interest in Kondo 
physic^i^E^El], jf the quantum dot is tuned in such a 
way that it carries a net spin, resonant tunneling leads 
to an incre ase of the conductance up to the unitary 
hmhPESEl. Por small dc-voltage bias V <^Tk the sys- 
tem can be described using linear response theory. How- 
ever, for V > Tk linear response theory starting from 
the equilibrium ground state is no longer applicable. 

In this paper we study a quantum dot in the Kondo 
regime (Kondo dot) with an applied magnetic field h in 
the regime max(y, T) ^ Tfc, where V is the dc- 
voltage bias and T the temperature. We diagonalize 
the Hamiltonian usin g infi nitesimal unitary transforma- 
tions (fiow equations p^^. Unlike in conventional scal- 
ing approaches the high energy states are not integrated 
out, instead the states are successively decoupled from 
large to small energy differences. Since all energy con- 



serving processes are retained, the steady current across 
the dot turns out to be included in the scaling picture. 
This current generates a decoherence rate F that cuts off 
the logarithmic divergences arising in the Kondo prob- 
lem, thereby making the situation max(y, \h\, T) 3> Tk a 
weak-coupling pr oblem. Previous ren ormalization group 
(RG) calculationj^ ^" l "7 | 28 | 29 | 3U | 3i | 3 2] ^j^g^^y estabhshed 

that decoherence effects due to spin relaxation processes 
play a key role in non-equilibrium. This was confirmed by 
a flow e quation analysis of the Kondo model with volt- 
age bia^^mMU^ Other new scaling approaches to non- 
equili brium pro blems like the real time renormaliz ation 
gj-Q^jJllIlIISlIlll g^jjj Coulomb gas representatiorP^BIl 

are consistent with this general picture and have added 
further insights. At this point one should also men- 
tion other new approaches like the scattering state nu- 
merical renormalization grou p^ , the time-dependent 
density renormalization grouppB the scattering state 
Bethe Ansat^^MS] and 1 /iV-expansion techniques^* that 
open up the possibility to describe the very challenging 
crossover regime for intermediate voltage bias V ^ Tx- 

In this paper we generalize the flow equation 
analysis'^'^ "^^ '^^ to include a magnetic field. A similar 
two-loop calculation based on the real time renormal- 
ization group was recently also performed in Ref."^. As 
main results we derive the spin-spin correlation function, 
the T-matrix and the magnetization in both equilibrium 
and non-equilibrium situations. Our results for the non- 
equilibrium static spin susceptibility xoiT, V) at zero ex- 
ternal magnetic field h — were already published in 
Ref.35 

Let us first have a closer look at the magnetization. 
The equilibri um magnetization is well known from the 
Bethe AnsatJ^. Previou s non-equilibrium perturbation 
theory calculation d^^ * '*^! -'^ for the magnetization derived 
the correct high voltage/temperature V,T 3> \h\ be- 
havior, but important logarithmic corrections in non- 
equilibrium are missing. Using the flow equation ap- 
proach up to two-loop order we will be able to calculate 
the magnetization including its leading logarithmic cor- 
rections consistently in the whole weak-coupling regime. 

The T-matrix and the closely related impurity s pectral 
function are also well studied objectsPEEMIMa^ 
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ertheless, some parameter regimes like combinations of 
magnetic field with nonzero voltage bias have not yet 
been investigated. We rederive the previous results and 
give additional insights into the crossover regimes. 

The equilibrium spin-spin c orrelation function is 
known in all parameter regimesj^^^^^^^nBll especially in 
the setting of the equilibrium spin boson model. We 
generalize our previous results'^* in both equilibrium 
and non-equilibrium to include nonzero magnetic fields. 
In addition we discuss in detail the interplay of the 
different dccoherence sources on the spin dynamics. 



The paper is organized as follows. In Sects. IB and IC 
we define the model and give a short introduction to the 
flow equation method. The flow equations for the Hamil- 
tonian and their scaling analysis are derived in Sect. [ll] 
(with additional details in appendices |A] and |b]). The 
transformation of the spin operator and the resulting cor- 
relation functions are shown in Sect. |III A| and Appendix 
[C] In Sect. IIIB we analytically derive the equilibrium 
zero temperature magnetization in leading order directly 
from the transformation of the spin operator. The cal- 
culation of the T-matrix is shown in Sect. IIII CI Numeri- 
cal results that show similarities and differences between 



voltage bias and temperature are discussed in Sect. HID 



B. Non-Equilibrium Kondo Model 

The Hamiltonian of a spin- 1/2 Kondo dot in a mag- 
netic field coupled to two leads is given by 



(1) 



h.C. 



where S is the impurity spin, a,f3 — l,r label the leads, 
a —'I, I is the spin index, ,. = ±V/2 is the chemical 
potential, and h is the magnetic field. Without loss of 
generality we assume V > 0. We are always interested 
in the isotropic Kondo model as is relevant in quantum 
dot physics, though most of our calculations can easily 
be generalized to the anisotropic case. 



Analogous to our previous calculationd^Hmni] split 
the operator space in even and odd combinations of 
fermionic operators from the left and right lead: 



fp<y 
9piy 



1 



1 



\/m? 

1 



1 



(2) 



VT+R 



where R = Ju/J^r- Note that the /- and g-operators 
obey fermionic anticommutation relations. If the Hamil- 



tonian ([T|) is derived from an underlying Anderson im- 
purity modeP^*^, the antisymmetric operators ffpo-'Spo- 
decouple completely from the dot and the Hamiltonian 
([ij can be written in terms of the /—operators only: 



H = J2^pflJpa-hS^ 



(3) 



P,(T 



((/pT/.T-/pV.i)^^ 

p.q 



def 



Ju + Jr 



where J 

JuJrr^^' The Hamiltonian 



and we have used Jf^ — J^i — 



3j looks formally like a 
standard Kondo impurity coupled to a conduction band, 
the only difference being the non-equilibrium occupation 
number distribution of the initial state derived from Q: 



nf{p) if la f pa) = 



, ep> V/2 



l+l/R 
1 



\ep\< V/2 (4) 
'V/2 



In equilibrium the Kondo temperature is given by Tk — 
Dy/pJ cxp(— where 2D is the bandwidth and p 
the conduction electron density of states (we assume a 
constant density of states). We will use this definition of 
the Kondo temperature in the remainder of this paper. 
For convenience we also set p = 1 in the following. By 
using the Hamiltonian ^ we will be able to describe the 
equilibrium and the non-equilibrium system in a unified 
scaling picture. For later reference let us already quote 
the result for the steady state current in the la rge d c- 
voltage limit for vanishing external magnetic fielcP^EIl 



T (T 



1 



V 



i?)(l + i?-i) Iu^V/Tk) 



(5) 



C. Flow Equations 



The flow equation methocP^^^ provides a framework 
to diagonalize a Hamiltonian using inflnitesimal unitary 
transformations. These are constructed using the differ- 
ential equation 



dH{B) 
dB 



^[il{B),H{B)] , 



(6) 



where the generator rj{B) is a suitable antihermi- 
tian operator. H{B = 0) is the initial Hamiltonian 
and H{B — oo) the diagonal one. The generic 
choice for the generator is given by the commutator 
r^{B) = [Ho{B),Hir,t{B)], where Ho{B) is the diagonal 
part of the Hamiltonian and 7?int(-B) the interaction 
part. With this definition of the generator one can 
define an energy scale Afeq = B^^^^ that corresponds to 
the remaining effective bandwidth: Interaction matrix 
elements with high energy transfer |Ae| > Afcq are 
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eliminated in the Hamiltonian ff(Afoq), while processes 
with smaller energy transfer are still retained. 

For generic many particle problems the flow generates 
new interactions, which appear in higher order of the 
interaction parameter. To keep track of the latter 
we introduce a parameter A = 1 in the Hamiltonian 
H = Hq + XHint- We will only take terms into account 
that couple back into the flow of the original Hamilto- 
nian up to a certain power of A. This corresponds to 
a loop expansion in renormalization theory: a n-loop 
calculation takes terms of order A""*"^ into account. We 
use normal ordering to expand operator products, see 
Appendix [A] for more details. 

To evaluate expectation values the operators have to 
be transformed into the diagonal {B = oo) basis. Any 
linear operator O is transformed using 

^ = [,(5), 0(B)]. (7) 



A generic operator will typically generate an inflnite 
number of higher order terms and one has to choose a 
suitable approximation scheme, which is again perturba- 
tive in the running coupling. 

The flow equation approach has been successfully ap- 
plied to various equilibrium many -body problems, like 
dissipative quantum system^21I13 tJig two-dimensional 
Hubbard model^^^, low-dimensional spin systems^" 
and stron g coupling models like the sine-Gordon 
modeP^I^ and the Kondo modeP^EH. It has also been 
successfully a pplied to n umerous non-equilibrium initial 
state problemsPMMIlSl. 



II. FLOW OF THE HAMILTONIAN 
A. Ansatz and Generator 



In the following we derive the Hamiltonian flow for the Kondo Hamiltonian (|3|. We use the ansatz 

Ho = J^epi/pUp. (8) 

p. a 

Hint ^ ■■ fl^h^ ■■ -4,{B) : ■)s' + \T. fl^hi ■■ S-+ : ^/.t : S+) 

+ E K,,rs{B){-- /Jt/^i^/.T : S-+h.c.) + E K!^,,rs{B){: flUfljsi : 5" -Hh.c.) 

p,q.r,s p,q^r,s 

where : : denotes normal ordering with respect to the system without Kondo impuritjl^, KpqlJ^{B = 0) = and 
Jpq^^^{B = 0) = J. For zero magnetic fleld the relations h ^ 0, jj^ = J^^ = J^^ — J^, and Kj^^ ^.^ = ^-^pg.rs — 
are fulfilled during the flow. The relations J^J^ — Jj^, Jj^^ — J^^, and if — ^srqp ^'^s always fulfilled 
due to hermiticity. An additionally generated potentia l scattering term is neglected since it has no influence on 
the universal low energy properties of the modeP^ISESl^ •\Ye also drop an uninteresting constant in the flow of the 
Hamiltonian. The straightforward derivation of the commutation relations yields the generator 

V{B) = \ E(ep - e,) (jl{B) : ^/.^ : -4,{B) : ^/.^ :) (9) 

+ \ T.(^P - + HB))4,{B) (: fl^Ul ■■ ■■ flifpl ■■ S+) 
p.q 

+ Y.('p-'i + 'r-^s + h{B))Kl^^,{B) (: fl^Ujl^fs^ ■■ S- - h.c.) 

p,q^r,s 

+ Y.('p~''i + 'r-^s + h{B))Kj,^^^,{B) (: fl^Ujljsi ■■ S- - h.c.) 

p,q^r,s 
p.q,r,s 

The resulting 2-loop flow equations are worked out in Appendix |B] The Hamiltonian is diagonalized in a controlled 
expansion if max(|/i|, V, T) 3> Tr-, which we assume in the following. Otherwise the running coupling becomes of 0(1) 
and an expansion in its powers is uncontrolled. 
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B. 1-loop Scaling Analysis 

The complete set of flow equations cannot be solved 
analytically due to the complicated momentum depen- 
dence. However, qualitative results for the low energy 
properties of the system can be worked out analytically. 
In the following we derive a simplified sca ling pict ure us- 
ing the so-called diagonal parametrizatiorP^^lESl 



Jpqi.B) - g-pq 



h{B)f 



(10) 



where pq — (cp + eg)/2. The energy diagonal equa- 
tions are easily obtained by setting Cq = Cp for the g^^^ 
terms and eq — ep + h for the g-^ terms. The diagonal 
parametrization can be seen as a generalization of the 
conventional IR-parametrization of scaling theory that 
allows for an additional dependence on the energy scale. 
Note that the flow of the running coupling does not de- 
pend on the momentum index p but on the energy scale 
ep-. We use Qp as shorthand notation for g^^. 

In the following we discuss qualitatively the flow of the 
1-loop equations. We find 



dh 
dB 



(11) 



iep-eq + h)ig^re~'^^^^-^^-^'^^' 



for the fiow of the magnetic field. Its small shift will be 
analyzed in the next section. 

For a nonzero magnetic field the Kondo couplings will 
not remain isotropic under the fiow. The running cou- 
pling for parallel scattering is given by 



dgl 

dB 



dgj 

dB 



J2{l-2nf{r)){ep-er + h) (12) 



x(.g^)2e-2S(^.-^^+'^)^ 



J2{l-2nf{r)){ep-er-h) (13) 

r 



and for spin-fiip scattering one finds 



dgp 

dB 



2 ^ 



(1 - 2nf{r)) ep 



^E(l-2n/(r)) I e„ 

r 



h 
~ 2 

h 

+ 2 

'2B{e 



(14) 



r + h/2f 



For convenience we generally drop the _B-argument of the 
running coupling and the magnetic field. At zero temper- 
ature the fiow of the running coupling can with very good 



accuracy be simplified using /(x) exp(— 2J5(x — c)^) ~ 
f{c) exp(— 2i?(x — c)^). This approximation removes the 
dependence of the running coupling and the summa- 
tions in ([T2|-( 14 1 can then be performed and lead to the 



following expression (p = 1 sets the energy scale) : 
de(l-2n/(e))(e + c)e-2S(e+c)^ ^ 

^ /g-2B(c-y/2)=^ ^-2B{c+V/2f\ 



2B \ l + R l + l/R 

This yields for parallel scattering 

dgl _ ( e-^B{e,+h+v/2r 



dB 2B \ l + R 

g-2S(ep+/i-V/2)^ \ 



1 + 1/i? j 

dgl _ (4-V2) fe-^Bie,-h+v/2f 
dB ~ 2B 1^ 1 + i? 

^2B(ep-/i-V/2)2 \ 



1 + 1/i? y ' 

and for spin-fiip scattering one finds 

dg^ _ 9p9i-h/2 ( e-2S(^P-V2+v/2)^ 



dB AB \ l + R 

g-2B(ep-h/2-y/2)2 \ 



1-l-l/i? 



5p<+V2 /e-2B(^P+'V2+V/2)^ 



4B \ 1 + i? 

^-2B(cp+h/2-V/2f \ 



l + l/R 



(15) 



(16) 



(17) 



(18) 



The fiow of the running coupling is cut off by an expo- 
nential decay unless Cp = —{h±V/2) for g^, ep = h + V/2 
for g^, or ep — +{h + V)/2 for gp . As a consequence 
the running coupling is strongly peaked at these energy 
scales. In the limit /i = this just corresponds to the 
strong-coupling behavior of the running coupling at 
the left and right Fermi level. The terms in 2-loop 
order cut off this strong-coupling behavior as we will 
see in the following sections. Replacing the exponentials 
in Eqs. (16l-(18l by 0-step- functions, these equations 



become equivalent to the p erturbative RG equations 
derived by Rosch et aZ.'SSHH. The different momentum 
dependence of the running coupling only leads to 
subleading corrections. 



5 



At nonzero temperature {T > 0,V = 0) one can un- 
fortunately not give a closed expression for 



de tanh 



(^)(e + c)e 



-2S(e 



(19) 



We therefore only discuss the asymptotic result for T ^ 
\h\. Since we are mainly interested in small energy scales 
£p ~* Oi we study the running coupling at the Fermi 
level only: g = gtJ2o^ ■ For B ^ T^^ the terms at 
high energies e S> T give the main contribution to the 
integral and we obtain the usual zero temperature scaling 
equatio 



dg 
dB 



2B 



(20) 



dg g' V2tt 1 



Note that Afcq = B~^/'^. For B > T"^ only energies 
e <C T contribute to the integral, since higher energies 
are cut off by the exponential. Therefore we linearize the 
tanh-function and obtain 



(21) 



This implies that the flow of the running coupling effec- 
tively stops for B > [tVb > 1). 

To obtain the numerical results shown later we have 
solved the 2-loop flow equations in diagonal parametriza- 
tion (10) since the solution of the full equations is very 



dB B 16 tVB 



resource intensive. We have verified in selected examples 
that this approximation agrees extremely well with the 
full set of equations. 



C. 2-loop Scaling Analysis I 

For the case of zero voltage bias and zero tempera- 

1(1) 

ture, the running couplings gp are strongly peaked 
at Cp = —{+)h, and g^ at ep — ±h/2. For a quali- 
tative analysis it is sufficient to replace the momentum 
dependent couplings by their peak values. For B ^ h~'^ 
we find the well-known 2-loop scaling equations for the 
anisotropic Kondo modePi' 



dg\\{B) _ gl{B) ,gi(B)5||(B) 



dB 
dg±iB) 



2B 

g±{B)gAB) 



4S 

g±iB){gm 



(22) 



gliB)) 



dB 2B 8B 

where g\\{B) = g\{B) = gi{B) and g^{B) - 
9±h/2^B)- The fiow parameter and the remaining ef- 
fective bandwidth are related by Afeq — B^^/^. The 
solution of Eqs. (22 1 is given by g{B) — gy±{B) — 

1 / hi{l / {y/BTx)) for an initially isotropic Kondo model. 
Additionally, we find a small shift of the magnetic field 



dh{B) 
dB 



§^V2^crf(V2BhiB) 



16B3/2 



(23) 



For small arguments (i? ^ /i 2) the error function is 
linear: 



dh{B) 
dB 



h{B)9l{B) 
2 2B 



(24) 



Using dg\\ldB = g\/{2B) Eq. ([24f yields: 
KB) 



ho exp 



2 \ln{l/{^TK)) HD/Tk) 



(25) 

where = h{B = 0). For large fiow parameters B ^ 
the error function is equivalent to the sign function. 
This yields a negligible additional shift of the magnetic 
field 



dhiB) _ gliB) 



2tt 



(26) 



dB 16B3/2 

of 0{ho/{ln{ho/TK))'^)- We can therefore use a constant 
magnetic field h{B) = h* for B ^ h^^, which is deter- 



mined from Eq. (|25| by h* ^ h{B = Hq"^). Notice that 



Bethe Ansatz calculation^l^l gnd a shift of the magnetic 
field to h* « /io(l — l/(21n(|/io|/7i<-)), which is consistent 
with our result (25 1 in the scaling limit D/Tk — > oo and 



confirms our approach. 

For B ^ we are then left with the fiow equations 
for the coupling constants 



dgwiB) 

dB 
dg±{B) 
dB 



= -gl{B)g\\{B)- 



/2T:\h* 



(27) 



= -9\{B) 



2Tr\h* 



16VB 



where we have neglected the 1-loop terms since they only 
contribute in 0{g^/B). The initial values are given by 



gyABo = ih*)-')=g* = 



1 

Hh*/TK) 



The differential equations ( 27 1 are solved by 

g 



g\\iB) 
g±iB) 



1 



-^wiVB- 

g* 



(28) 



(29) 



1 + T^i^-^) 



with 



Fii =ri 



/2n 



{g*r\h* 



(30) 



In Sect. Ill A we will see that up to a pref actor Fy can be 
identified with the longitudinal and Fj^ with the trans- 
verse spin relaxation rate: 



1 



OC Fii 



1 



OC F I 



(31) 



We will now already take this identification for granted 
so that we can compare our result ( 30 1 with literature 
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values. The spin relaxation rate in the limit that the 
thermal energy is much smaller than the magnetic en- 
ergy was first calculated iiP2l using unrenormalized per- 
turbation theory: it agrees with (30) if one uses the same 
approximation g* — > g. Ref.^^ also derived Ti = T2/2 in 
this limit, which in our calculation is hidden in the obser- 
vation that the proportionality factors in ( |31[ ) differ due 
to the different decay laws in ( 29 1 . We will not analyze 
this in more detail here since we will later in Sect. Ill 
even calculate the full line shape of the dynamical spin 
susceptibility. 



D. 2-loop Scaling Analysis II 

So far we could use the peaks of the running coupling to 
derive a simple scaling picture in equilibrium. Applying a 
dc-voltage bias yields a splitting of these peaks by ±V^/2 
since the resonances are pinned to the Fermi levels. As 
shown in our previous 

calculatioiJ221M135] ^^^s can be taken 
into account on the rhs of the flow equations by averaging 
over the splitting of the peaks, e.g. 



h+V/2 



V 



de gi{B) 



h-V/2 



In the previous section we expanded the flow equations 
for small flow parameter B <^ h^"^ and for large flow 
parameter B ^ h~'^ . If a dc-voltage bias is applied we 
flnd four energy scales that determine small and large 
B, namely |y -I- \V — h\, V and \h\. So in principle 
we would have to discuss the flow equations separately 
in all five regimes of the fiow. However, we can restrict 
the following discussion to the initial flow and the flow 
at very large flow parameter B ^ B^ — Aq^, where 
Aq = min {|V^ -I- |F — h\, V, \h\}. One can numerically 
verify that the flow in the intermediate regimes only leads 
to small corrections. 

In equilibrium at nonzero temperature we can 
again use the peaks of the running coupling to an- 
alyze the flow. Here the relevant energy scales are 
given by T and h coth{h / {2T)) and we deflne Aq = 
min {T,/i coth(V(2r))}. 

For small flow parameter B (initial flow) we find the 
usual scaling equations ( 22 1 and a small shift of the mag- 
netic field. In the regime B ^ Aq ^ the 1-loop terms are 
negligible. We are then left with the fiow equations 



dg^iB) _ 9l{B)g±{B) gl^B) 



dB 

dg\\{B) 
dB 



2Vb 
9i{B)gi\{B) 



B 



ci 



C2 



2Vb 



C2 



(33) 



The initial values are .g||/_L(-Bo) = 9\\/± the constants 



are given by 
Ci{h*,V) = 
C2ih*,V) = 



/2^ V 

IT {1 + R){l + 1/R) 

/2^\V + h*\ + \V-h* 



(34) 

\h*\{R+l/R) 



{l + R){l+l/R) 



in non-equilibrium (V ^ 0, but zero temperature T = 0). 
In equilibrium (V = 0, T ^ 0) we have 



ci(/i*,r) 

C2{h*,T) 



/27r, 



4 
27r 



-T 



(35) 



h* coth 



2T 



where we have defined h* = h{Bo). One easily shows 
that the solution of 



dg\\{B) 
dB 



9±{B) 



gl{B) 9\{B) hg^^Y ^ \ 

-ci - I ^^C2 - .giici (36) 



B 



B 



\ 



gl{B)^+g\AB) 



{9lY 



C2 



C2 



(32) also solves Eqs. (ISSl). Thus Eqs. {Isel) are an equiva- 



lent formulation of Eqs. (33 1. The remaining fiow equa- 



tion is an Abel differential equation of the first kind 
whose general analytic solution is impossible. There- 
fore only asymptotic results can be obtained. Note that 
the coupling constants asymptotically flow to zero or to 
a nontrivial flxed point g\\{B) = gJJ — C2{g*i)'^ / {g*<^ci) 
in Eq. ( 36 ) . This nontrivial flxed point can only be 



reached in the anisotropic Kondo model with initial val- 
ues g^{B = 0) < \g{{B = 0)| and \h\ < V,T. It is 
an unphysical artefact of our approximations in the flow 
equation calculation, where higher order terms need to 
be included in this regime. We will not say something 
anything about this part of the parameter space of the 
anisotropic Kondo model for the remainder of this paper 
and now return to the isotropic model. 

The equilibrium zero temperature behavior has already 
been analyzed in the previous section. Notice that the re- 



sult ( 30 1 for the spin relaxation rate derived there holds 



generally when the magnetic energy dominates, that is 
for large magnetic field h ^ V,T. We next look at 
the opposite limit of vanishing magnetic field, h — 0. 
The Kondo model remains isotropic and longitudinal and 
transverse relaxation rates coincide. Eq. (36 1 is solved by 



9iB) = 



1 



.9' 

t{Vb 



(37) 



where g{B) — g\\{B) = g±{B) with the spin relax- 
ation rate F — A{g*)'^ci. The equilibrium finite tem- 
perature spin relaxation rate therefore shows the ex- 
pected Korringa-behavior proportional to temperaturd^. 
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while the non-equilibrium zero temperature spin relax- 
ation rate is proportional to the voltage bias (or more 
accurately according to (34) and (jsj: proportional to the 
current across the dot), which agrees with the perturba- 
tive non-equilibrium RG-result^''. 

For intermediate values of the magnetic field we find a 
competition between the quadratic and the cubic terms in 
the running coupling in Eq. ( 36 1 . For very small energies 
(_B — !■ cx)) the quadratic term dominates the flow if it 
is nonzero and the running couplings decay like g\\ ~ 
B^^/^ and g_L ~ If the cubic term dominates 

both running couplings are proportional to B~^/^. One 
can still analyze this analytically for small magnetic field 
y, T 2> \h\: Eqs. (33) are approximately solved by 



5||a(S) 



(38) 



i + r||/^(VB- VSo) 



since « g*^_ and ci « C2. In the high voltage regime the 
energy scales where the algebraic decay of the couplings 
sets in are given by 



MglY 



V 



(i + i?)(i + i?-i) 



{gl? 



W- 



+ \V ^h*\ + \h*\{R + R-^) 

(l + i?)(l + i?-l) 



Here we have kept the leading /i-dcpendence to show that 
as expected only the spin flip coupling sees both the mag- 
netic fleld and the voltage bias in its relaxation rate. We 
flnd similar behavior in equilibrium at nonzero tempera- 
ture for < T: 



^\\{h\T) = 
T^{h\T) = 



(40) 



{g±) h* coth 



h* 

2r 



We have now analytically derived the qualitative scal- 
ing behavior and the related energy scales of the Kondo 
model as a function of temperature, voltage bias and 
magnetic fleld. In particular we have seen how sufficiently 
large temperature, voltage bias (more accurately: a suffi- 
ciently large current ([S])) or magnetic field can make the 
Kondo model a weak-coupling problem, where the cou- 
pling constants decay to zero and therefore allow for a 
controlled solution using flow equations. Similar observa- 
tions have be en m ade using other renormalization group 
techniques'^ 38llo]^ ^j-^g next chapter we will turn to 
a completely numerical solution of the flow equations in 
order to obtain quantitative results. This is also the only 
way to analyze the behavior when the magnetic fleld is 
of the same order as the voltage bias or the temperature, 
which was excluded in the above analytical discussion. 
Still, the analytical results obtained so far are important 
because they will serve as our guidelines to understand 
and interpret the results of the numerical solution. 



III. RESULTS AND DISCUSSION 

We restrict the following discussion of numerical re- 
sults to symmetric coupling to the leads, R = 1. The 
extension to i? 1 is straightforward and corresponding 
results can be obtained without further complications. 
We also want to mention that the approximations in our 
calculation in this paper do not allow us to obtain a more 
accurate result for the current than the previously known 
expression (|5| , and therefore we will not elaborate on the 
evaluation of the current in this chapter. 



A. Spin-Spin Correlation Function 

Since the interacting ground state becomes trivial in 
the B = oo basis, we transform all operators into the 
diagonal basis before calculating their expectation values. 
We make the following ansatz for S^: 



(39) S'iB) = h%B)S' 



M{B) 



(41) 



'Y.^,,{B){:fl^U-S-+:flUr-S^), 



where h^{B = 0) = 1, -fpq{B = 0) = and M{B = 
0) = 0. For the transformation of the spin operator it is 
sufficient to use only the first order part of the genera- 
tor ([9|, that is to neglect terms in O(J^) in the genera- 
tor: In Ref.lSSIwe showed that this approximation already 
yields results including their full leading logarithmic cor- 
rections. 

The decay of the spin operator into a different struc- 
ture under the unitary flow is described by the flow equa- 
tion for the coefficient h^{B): 



= -Y.(^Ap) + Ml)-^Mp)nA<l)) (42) 
x{ep-e, + h)J^^{Byf,,{B) . 
For the flow of the newly generated c-number we flnd 



dM{B) 
dB 



= J2inf{p)~nf{q)){ep-e, + h) (43) 



p,q 



For zero magnetic fleld the relations J^^^B) — J^(i?) 
and jpq{B) = —jqp{B) arc fulflUed. Using these relations 
one easily shows M{B) = 0. We will later see that M is 
just the magnetization and therefore it makes sense that 
M only becomes nonzero during the flow if an external 
magnetic fleld is applied. 

The flow of the newly generated operator structure in the 
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spin operator is given by 



dB 



(44) 



r 

x-frg{B) + {er~e,)4^{Bhpr{B)) . 

In the sequel we will focus on the longitudinal spin sus- 
ceptibility. The calculation of the transverse part follows 
exactly the same route and the transformation laws of 
are given in Appendix [c| 
During the initial flow h^{B) is nearly unchanged: its 
flow is only of 0(1711(5 = 0) — gi\{B)). For simplicity we 
restrict the following discussion to equilibrium and zero 
temperature; the extension to V, T > is straightfor- 
ward. In lowest order the solution of Eq. (44 1 is formally 
given by 



-h) JdBi h^Bi)J^{p,q,Bi) 



(45) 



With Eq. (42 1 follows (using diagonal parametrization) : 



1b 



(46) 



where h is the external magnetic fleld. Notice the simi- 
larity to the 2-loop flow equation for t^y (27 1. Therefore 
the 5'^-operator begins to decay on the same energy scale 
at which the strong coupling divergence of g\\ is cut off. 
By a similar argument one can show that the decay of 
the S'^/^'-operators is related to the flow of gj^ : the decay 
starts on the same energy scale that cuts off the strong 
coupling divergence of g±. We conclude that the energy 



scales Fy and F_l determine the decay of the spin opera- 
tors parallel and perpendicular to the external magnetic 
field. 

It is this observation that relates the decoherence rates 
Fy and r_L to the physical spin relaxation rates: 1/Ti and 
1 /T2 are defined through the broadening of the resonance 
poles in the longitudinal and the transverse dynamical 
spin susceptibilities'^" . Now for B ~ ^^^± ^lU excitations 
with energy transfer much larger than the decoherence 
rate are integrated out. Since the spin operator has not 
yet decayed on this B-scale, the broadening of the res- 
onance pole can therefore not be larger than the corre- 
sponding energy scale. The algebraic decay of the spin 
operator just corresponds to the broadening of the res- 
onance pole. Hence (up to a prefactor) l/Ti ~ F|| and 
I/T2 ~ Fj_. The width of the broadening of the reso- 
nance poles is therefore (up to an uninteresting prefac- 
tor) automatically given by the decoherence rates defined 
in Sects. |IIC| and |IID[ As already discussed there, our 
results for the longitudinal and the transverse spin relax- 
ation rates agree with previous results in the literature 
in their appropriate limits. 

The symmetrized spin-spin correlation function is 
defined as C"(ii,i2) = i({S'°(ii), ^"(ts)}) and 
the response function as Xa(ii:i2) = —iQ(ti — 
t2){[S'^{ti), S'^{t2)]) ■ In equihbrium these expectation 
values are evaluated with respect to the ground state 
or thermal state defined by H{B — 00).'^ In non- 
equilibrium they have to be evaluated with respect to 
the current-carrying steady state. This was discussed in 
Ref.'^^: By explicitly following the time evolution of the 
initial state until the steady state builds up, we could 
show that we can work with the same state as in equi- 
librivim in the present order of our calculation (notice 
that this does not hold for the evaluation of the current 
operator itself). 



The Fourier transformed S^-S^ correlation function is therefore both in equilibrium and non-equilibrium given by 



r^^^^ E - M^P + ^ + h)) + l'^nnf{e,){l - n,{e, + h))) (47) 



7r(l -t- sgn(/i) 



where the tilde denotes the value at i? = 00. The corresponding imaginary part of the Fourier transformed response 
function is 



7r(l - sgn(/i)) 



' ep,ep—uj-\-h 



rnf{ep){l-nf{ep^u; + h))) (48) 



p 



7r(l + sgn(/i) 



p 

I 



the real part is accessible via a Kramers-Kronig transfor- 



mation. The correlation function is a symmetric func- 



9 



tion of w, the imaginary part of the response func- 
tion is antisymmetric. Both functions do not depend 
on the sign of h. In equihbrium the fluctuation dissi- 
pation theorem-^ relates the imaginary part of the re- 
sponse function and the spin-spin correlation function 
by Xzi^) = tanh(ijj/(2T))C^(a;). In non-equilibrium the 
fluctuation dissipation theorem is violated in general. For 
completeness the corresponding expressions for S'^/?' are 
given in Appendix [Cl 

Typical equilibrium zero temperature spin-spin corre- 
lation functions are shown in Fig. IT] At zero frequency 
we find a 5-peak with strength M {B = cx))7r/2 in the 
correlation function (47) due to the nonzero spin expec- 
tation value (it is not plotted for obvious reasons). For 
convenience we assume ft, > in the following discussion. 
The maximum of the spin-spin correlation function (ig- 
noring the (5-peak at w = 0) is as expected at a; ~ /i* and 
it decays with increasing magnetic field (see the inset of 
Fig.[T]). For jwl < ft* the correlation function vanishes ex- 
actly in the present order of the calculation, for \u!\ ^ |ft| 
we find C^(w) ~ l/(|w|(ln(|w|/Tjf))2). Fig. |2] shows 
the buildup of this characteristic behavior of the spin- 
spin correlation function also for nonzero voltage biad^ 
upon increasing the magnetic field. The inset shows the 
corresponding plot in equilibrium for nonzero tempera- 
ture. In non-equilibrium we find pronounced peaks at 
\uj\ ~ \h-V\, |w| - \h+V\, and - ft for ft > l^. The 
peaks at w ~ ±|ft — y| join for h < V and build up the 
zero frequency peak. 




FIG. 1: Equilibrium S^-S^ correlation function for various 
magnetic fields, V=T=0. The inset shows the decay of the 
peak height, which approximately follows a power law. 










— h = 


20 




0.006 




.... h = 


22 \ 








-- h = 


24 Tk 








.-. h = 


26 


H 






— h = 


28 




0.004 














.... h = 


30 














0.002 









Notice that for nonzero temperature all additional 
peaks are smeared out (inset of Fig. [2]). This exemphfies 
a key difference between non-equilibrium and nonzero 
temperature that will keep reappearing in other dy- 
namical quantities. The non-equilibrium Fermi function 
Q retains its characteristic discontinuities, which 
lead to strong-coupling behavior yielding Kondo-split 
peaks in dynamical quantities. These peaks are only 
cut off by the decoherence rate and not by voltage or 
temperature itself, and therefore much more pronounced. 



FIG. 2: Emergence of the large magnetic field behavior of the 
spin-spin correlation function depicted in Fig. [l] for a non- 
equilibrium situation {V = 20Tk,T — 0). The inset shows 
typical equilibrium correlation functions for nonzero temper- 
ature T. Here all features except the zero frequency peak are 
smeared out. 



B. Magnetization 



The spin-spin correlation function of the Kondo Model 
in a magnetic field has so far mainly been studied in the 
context of the spin boson model^" For high magnetic 
fields no previous results exist, since high frequencies 
are difficult to access by numerical methods like NRG. 
Paaske et al^ studied the transverse dynamical spin 
susceptibility for high voltage bias, which can be calcu- 
lated within our approach using the transformation of 
the spin operators perpendicular to the magnetic field 
in Appendix [Cl Using a Majorana fermion representa- 
tion Mao et alA^ obtained the low frequency properties 
of this correlation function in the case of dc- voltage bias 
and nonzero temperature in agreement with our results. 



From the ansatz (41 1 follows directly that the magneti- 
zation of the dot spin is given by M(B = oo). However, it 
turns out that the operator decays slowly with B and 
therefore also the magnetization converges slowly, mak- 
ing the analysis difficult. Still, we can use the following 
trick to obtain the leading behavior of the magnetization 
for T = V = analytically. Clearly 



2{S') = 2ft^(c5o)(c5o|S'^|oo) 4- Af(oo) 
= ft^(oo)sgn(ft(oo)) -I- Af(oo) , 



(49) 



where \oo) is the ground state of Hq{B = oo). Note that 
ft^(oo) = 0. For convenience we assume ft > in the 
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following. We rewrite Eq. (49 1 to the form 

d{h'{B) + M{B)) 



2{S'') = /i^(0) + Af(0)+ / dB 



1-2 fdB J2nf{q){l 

n P'l 



V 



dB 



x{ep-eq + h)Jpg{B)jpq{B) 
Using the parametrization 



9±iB) 



2(ep - 



h) 



we find 



2(5'^ 



^5 



(50) 



(51) 



(52) 



where f{B) = 1 for B < h'"^ and /(B) = for B > /i^^. 
Neglecting higher order corrections we find 



— analytic result 
o numerics 




h= 100 T 



D/T„ 



FIG. 3: The magnetization as a function of the external mag- 
netic field calculated from the analytic result ( 54 1 and from 



the numerical solution of the flow equations for D — IQ'^Tk- 
The inset shows the bandwidth dependence of the magneti- 
zation at ft = lOOTff. 



nonzercPni xhc imaginary part of the T-matrix is given 



Im[T<,(c^)] 



dte(t)({O,(t),Ot(0)})e-* , (56) 



2{S'' 



1 - 



dB 



4B 



(53) 



With dg\\{B)/dB = g']_{B)/{2B) follows 



2(5" 



1 - 



21n(VrA') 21n(D/Tx) 



(54) 



which (for D/Tk ^ oo) is to leading order the Bethe 
Ansatz resulli^. Fig. [s] shows the good agreement be- 
tween the analytic expression and numerical results for 
high magnetic fields. For fields of 0(10 T^) we see de- 
viations from the analytic result due to the perturbative 
nature of our approach. The inset shows the bandwidth 
dependence of the magnetization in good agreement with 
Eq. pll). 



Unfortunately, a similar analytical calculation for 
T^V > Q has not been possible. We will present nu- 
merical results in Sect. IIIID 31 



where 

0^{B) = Y.^Ut{B)h^S- + Ul{B)h^S^) (57) 
fe 

0^{B) = Y.{V^^{B)h^S+ -Vl{B)h^S^) , 

k 

U^'\b = 0) = J-L/T/2 and V^'^{B = 0) = In 
lowest order the flow equations for the spin up component 
are given by 



dB 



J ^(1 - 2n/(r))(ep - + h)U^J^^ , (58) 



dU^ 
dB 



4 ^ 

r 

4 ^ 



{l-2nf{T)){ep-er + h)Uljtp 
{\-2niir)){ep~er)U^Jip. (59) 



C. T-Matrix 

The scattering of conduction band electrons from lead 
a to lead /3 is described by the T-matrix Tap^ai'^)- It is 
defined via the electron Greens function 

If the Hamiltonian ([l]) is derived from an Anderson im- 
purity model only one eigenvalue of the T-matrix is 



Comparing the latter equations with Eqs. (12 1 and (14 1 



one already notices their similarity to the fiow of the 
running coupling in 1-loop order. In the following we 
work out the details. Using the approximations from 
Sect. El] we find 

dUl 9i+h,2K+h /e-^(^^+''+^/2)' 



dB 



2B \ 1 + R 

g-s(£p+h-y/2)^\ 



1 + 1/i? 



(60) 
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FIG. 4: Imaginary part of the spin-down T-matrix as a func- 
tion of the magnetic field (V = r = 0,D = 10^r;f). The inset 
shows a comparison with the Bethe Ansatz result for uj = 0. 



Neglecting a factor two i n t he exponential, this equa- 
tion is equivalent to Eq. (16 1 provided Uj, = gj,/2 and 



+h/2 = 9p 1'^- Analyzing the flow of the spin-flip com- 
ponent we find 



-h/2 



-B{ep-h/2+V/2y 



dB 



p-Biep-h/2-V/2f \ 



l + l/R J 

l + R 

g-B(ep+h/2-V/2)2 \ 



■9ep + fe/2'^ep + fe/2 I 6 

Ib I " 



l + l/R 



(61) 



Again neglecting the factor two in the exponential, this 
equation is equivalent to the 1-loop flow equation ( 18 ) 
for . One easily shows that higher order terms in the 



transformation of the operator ( 57 1 have the same effect 



on the flow as the 2-loop terms in the transformation of 
the Hamiltonian. The calculation for nonzero tempera- 
ture is again more difficult, nevertheless we find the same 
relations between the flow of the operator and the run- 
ning coupling. 

Doing an analogous argument for the y-terms we iden- 
tify 



Ul{B) 
V^{B) 



am 



v,^{B) 



-h/2 



(62) 



-h/2 



(B) 



Therefore the operators completely decay into more 
complicated objects for B oo. Since calculating the 
latter is resource intensive (three momentum indices), it 
is more economic to evaluate the T-matrix at the de- 
coherence scal^, where the decay of the couplings sets 




FIG. 5: Decay of the zero frequency peak in the spin-spin 
correlation function with increasing magnetic field, a) = 
20Tk, b) T = 20Tk. 



in and higher order terms in the transformation of the 
observable are not yet important: 
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(63) 



Here the hat denotes functions at the decoherence scale. 
Though the further flow leads to a decay of O^, the spec- 
tral function remains unchanged for B > rj7?, , where 



Th/± is the dominant decoherence scale. In (63 1 we can 



replace the expectation value of at the decoherence 
scale by the magnetization of the system since the 
operator decays noticeably only for B ^ '^\\^±- 

As suggested by Rosch et al^ we use Fermi functions 
broadened by the decoherence scale rj_ to describe the 
distribution function for the /-operators at the decoher- 
ence scale nf{uj). This avoids the costly full numerical 
solution to i? ^ CX3 and yields results that are virtually 
identical. In equilibrium at small temperature T ^ \h\ 
the distribution function is then given by fif(uj) = frii^), 
where /r(w) = 1/2 — arctan(ti;/rj_)/7r. At high temper- 
ature T ^ \h\ the spin expectation value (S^) vanishes. 
Then the distribution function only enters in subleading 
order. Note that the imaginary part of the T-matrix in 
general depends only weakly on the details of the broad- 
ening scheme. In non-equilibrium the step functions at 
both chemical potentials have to be broadened yielding 
hfiij) = fr{oj + V/2)/{l + R) + fr{i^-V/2)/{l-^l/R) for 
the distribution function. The additional factor of 1/4 in 
comparison with the result obtained by Rosch et alJ is 
due to our different deflnition of J. For symmetric cou- 
pling R — 1 the spin-up and the spin-down component 
are related by Im[T|(tj)] = Im[T|(-tj)]. 

The imaginary part of the T-matrix and the spectral 
function are related by A^(lo) = —lin[Ta-{uj + id)]/7r. 
Fig. |4] shows spectral functions for several values of the 
magnetic fleld. They are strongly peaked at w ~ /i. 
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FIG. 6: a) Kondo splitting of the sharp edge in the spin-spin 
correlation function (magnetic field h = IOOTk) due to a small 
voltage bias, b) Small temperature only leads to broadening 
and no Kondo splitting. 



Rosch tt al^ studied this structure in detail by analyzing 
the spectral function normalized to 1 as a function oiuj/h. 
Since we have included the shift of the magnetic field, we 
will do likewise as a function of . In agreement with 
the results derived by Rosch et aZ.'^lwe find that the width 
of the left flank is approximately proportional to the de- 
coherence rate (30 1, leading to a sharpening of the left 



flank for increasing h, while the width of the right flank 
increases for increasing h. 

The imaginary part of the T-matrix at zero fre- 
quency is relat ed t o the magnetization via the 
Friedel sum rule.'^^'^^ Inserting the leading term of 
the Bethe Ansatz result'SS one finds Im[T;,(0)] = 
-sin^ (7r/(41n(VTV))) /vr. The inset of Fig. g shows a 
comparison between the Bethe Ansatz and the flow equa- 
tion result. Again we find very good agreement for high 
magnetic fields and deviations for fields of 0(10 Tk)- 

For large frequencies the spectral function decays pro- 
portional to 1 / {ln{uj /Tk))'^ J ^ ^ which is consistent with 
our results. Also Bethe Ansatz calculation^^^ show that 
the maximum of the spin-down spectral function is at 
h* « h{l — l/(21n(/i/7if)), which is consistent with 



our shift of the magnetic field ( 25 1 in the scaling limit 
D/Tk ~> oo. 



D. Voltage Bias vs. Temperature 



1. Spin- Spin Correlation Function 



The spin-spin correlation function at zero ma gnetic 
field {V OT T ^ Tk) shows a zero frequency pealP^*^. 
In Fig. |5] we show its decay due to an applied magnetic 
field. Again the zero frequency S-peak in Eq. (47 1 is not 



FIG. 7: Sum of both spin components of the T-matrix at 
zero magnetic field for various values of a) voltage bias and 
b) temperature. 



plotted. The sum rule 



dLoC'iuj)^ M^-+ I dioC^iuj) (64) 



is not fulfilled exactly since we neglect higher order terms 
in the transformation of S^. The error is typically of 
order one percent. Here C^{lu) denotes the 7pg terms in 
Eq. (liTl). 



For increasing magnetic field the magnetization M in- 
creases. Due to the sum rule and the fact that C^(w) is 

a non-negative function, an increase of M must lead to a 
decrease of C^(a;), leading to a decay of the correlation 
function for lu 0. 

At first glance the decay of the zero frequency peak 
looks similar for both the equilibrium and the non- 
equilibrium case. Only the relative decay of the maxi- 
mum as a function of V/h and T/h seems to be different. 
On closer inspection we find additional peaks at ^ \h\ 
for V > \h\: their height increases with the magnetic 
field, see Fig. [2j In equilibrium for nonzero temperature 
these peaks are smeared out. For high frequencies we 
find the usual C^(w) ~ l/{\uj\{ln{\uj\/TK))^) behavior. 

In Fig. [6] a) we show the Kondo splitting of the sharp 
edge at Lu = h* in the correlation function due to an 
applied small voltage bias. The two new peaks are lo- 
cated at \uj\ ~ \h* ± V\. On the other hand, for small 
temperature we again only find a broadening effect. 



2. T-Matrix 

Fig. [7] depicts the sum of both spin components of the 
T-matrix for vanishing external magnetic field, that is 
the full impurity spectral function. For nonzero voltage 
bias {V ^ Tjf) and zero magnetic field one finds the 
characteristic Kondo split peaks at w ~ ±F/2. For 
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FIG. 8: a) Magnetic field shift of the Kondo split zero fre- 
quency peak of the spectral function for V = 20Tk, T = 0. b) 
Nonzero temperature in equilibrium (here: T = 20Tk, V — 0) 
leads to broadening. 




FIG. 9: a) Splitting of the peak in the spectral function due 
to a small voltage bias (external magnetic field h = IOOTr-). 
b) Again temperature only leads to broadening. 



nonzero temperature one only observes the expected 
broadening of the zero frequency peak. These observa- 
tions are consistent with the results obtained by NRGP 
and perturbative RG*^. Applying a small magnetic 
field leads to a shift with the magnetic field strength 
and an asymmetric deformation of the peaks. Typical 
curves are shown in Fig. [8] For large magnetic fields we 
have already discussed in Sect. Ill C how the spectral 
function of the equilibrium zero temperature Kondo 
model develops a pronounced peak at oj ~ /i, see Fig. |4] 
In Fig. |9] a) we show the splitting of this peak into two 
peaks at a; ~ h±V/2 due to a small voltage bias. Again, 
applying a small temperature only leads to a broadening 
of the peak. 
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FIG. 10: The magnetization as a function of the external 
magnetic field for a) fixed temperature and b) fixed voltage 
bias for D = IOOTk- For increasing temperature one ap- 
proaches the free spin behavior as expected. For increasing 
voltage bias one approaches the non-equilibrium perturbation 
theory result based on a quantum Boltzmann equation plus 
rate equations (QB/REJ^, which becomes exact in the limit 
V/Tk oo. 



difficult using NRG. For example the finite temperature 
broadening in Fig.[6]and our results for the T-matrix with 
magnetic field plus voltage bias or nonzero temperature 
have not been previously obtained using other methods. 



3. Magnetization 



In principle the magnetization can be extracted from 



We can see that it is straightforward to resolve sharp 
features in the dynamical quantities at large frequencies 
using the flow equation approach, which is notoriously 



the spin-spin correlation function via the sum rule ( 64 1 . 
However, due to the approximations in our calculation 
the sum rule is not exactly fulfilled and we were only 
able to extract qualitative results via this route. More 
accurate results can be obtained by analyzing the flow 
of M{B) directly. We were able to reproduce previously 
known results from Bethe Ansatz and non-equilibrium 
perturbation theory. 

In equilibrium the exact magnetiz ation is accessible 
by solving the Bethe Ansatz equation Assuming 
h > the asymptotic results relevant for this paper are 
given by the zero temperature magnetization M(h,T — 
0) = 1 - l/{21n{h/TK)) for h > Tk, and the high 
temperature magnetization M{h,T) — tanh(/i/(2r)) for 
T 3> Tx and T ^ h. The high temperature result is of 
course just the magnetization of a free spin. 

Previous non- equilibrium perturbation theory 
calculationd^lElIlll jn the limit F » T/^ or > 
found Mpt{h,V) = 4:h/{2\h\ + \h + V\ + \h ~ V\) for 
the magnetization. Here the important logarithmic 
corrections at zero voltage bias are missing since 
A/pt(/i,l/) =sgn(/i) ioTV<\h\. 

In Sect. [TTTB] wc have already derived the zero temper- 
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FIG. 11: Crossover from the equilibrium zero temperature 
magnetization a) to the high temperature and b) to the high 
voltage result for D = IOOTk. 



Refs.EIHl). 

As a final comment we want to mention that while the 
flow equation approach has to rely on numerical evalua- 
tions of complicated sets of differential equations (at least 
if one is interested in quantitative results beyond leading 
order) , it does allow one to study all combinations of the 
parameters voltage bias, temperature and magnetic field 
in one framework. As an outlook this should be useful 
for investigating other more complex quantum dot struc- 
tures in the future. 
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ature magnetization within the flow equation framework. 
Figs. [To] and [TT] show the crossover between the equilib- 
rium zero temperature result and the asymptotic high 
temperature result or the asymptotic high voltage bias 
result. These crossovers are smooth and show the ex- 
pected reduction of the magnetization for finite < 
that are missing in previous calculations. It should be 
noted that there is a noticeable dependence of the results 
in Figs . [T0| and [TT] on the bandwidth similar to Fig. [3] 

IV. CONCLUSION 

In this paper we have employed the fiow equation ap- 
proach to derive a consistent scaling picture of the equi- 
librium and non-equilibrium Kondo model in its weak- 
coupling regime. The weak-coupling regime is realized for 
sufficiently large voltage bias V, magnetic field h or tem- 
perature T as compared to the equilibrium Kondo tem- 
perature: max(/, \h\, T) ^ Tr-, where / is the current ([5| 
across the dot. Our calculation allowed for the evaluation 
of static and dynamical quantities including their leading 
logarithmic corrections. Specifically, we have studied the 
spin-spin correlation function, the magnetization and the 
T-matrix as functions of V, h and T and explored various 
crossover regimes. 

As emphasized by Millis et al) ^" * "^^ !, the non-equilibrium 
noise generated by the steady state current across a 
quantum impurity can to leading order be approximated 
by thermal noise with an effective temperature roff(-^), 
but with important differences between non-equilibrium 
noise and thermal noise occuring beyond leading order. 
We could see this explicitly in many of our dynami- 
cal quantities, where non-equilibrium conditions due to 
a voltage bias lead to effects like Kondo splitting and 
strongly enhanced logarithmic corrections (this is also 
especially noticeable in the static spin susceptibility, see 



APPENDIX A: NORMAL ORDERING 

In this section we briefly sum up some properties of 
normal ordered operators that are frequently used in flow 
equation calculations. For more details we refer to ReiW 

In the following Ap denotes creation and annihilation 
operators, the a's are c-numbers and P{{Ap\) is a prod- 
uct of operators from the set {^p}. The rules for Wick's 
normal-ordering are given by: 

1. Numbers are unchanged: 

: a:= a (Al) 

2. Normal-ordering is linear 

: aiPi({Ap}) + a2P2({Ap}) : = (A2) 
ai : Pi({Ap}) : +a2 : P2{{Ap}) : 

3. Recurrence relation 

A, : P{{Ap}) : = : A,P{{Ap}) : (A3) 
^Y-^ . dP{{Ap}) ^ 



where 

C,r = mA^Ar\^) (A4) 

for a pure reference state \^) or 

Cqr = Tr(pA,A,) (A5) 

for some mixed state described by the density matrix 
p. Typically the ground state of the non-interacting 
Hamiltonian is chosen as reference state \^) . 
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From the recurrence relation (A3) one can derive 
Wick's first theorem 



A A 



(A6) 



In the 2-loop calculation we neglect terms with four or 
six fermionic operators on the rhs. in the following since 
these terms would enter the calculation only in 3-loop 
order. We again need the commutator 



^P„_l ^ ] Cpn-iqn-l 



An 



From this relation follows that the commutation of neigh- 
boring fermionic operators picks up a minus sign, bosonic 
operators commute. The product of two normal-ordered 
objects can be calculated from Wick's second theorem. 
The fermionic version is given by 



c\,ci :, : 4,C24,C3 :] = 

5i',2'5i,2'Wl')-n(l)) :4'C3 
-<^i',2<5i,3'Wl')-«(l)) 



-<5i',3'5i.2'Wl')-"(l)) :4'C2 
+(5i',3'5i,3'("(l') -"(!)) : 4'C2 



(B4) 



: P,{{A,}) :: P^HA,}) : = : exp (E ^"-^^^ )(A7) 



-.P,{{Ap})P2{{Bp}) 



A=B 



APPENDIX B: TRANSFORMATION OF THE 
HAMILTONIAN 

The derivation of the flow equations for the Hamilto- 
nian (|8| is straightforward. Only some preliminary rela- 
tions are needed. Products of spin operators are easily 
calculated using the standard spin operator algebra. The 
relations 



[AS-,BS''] 
[AS+,BS'] 
[AS+,BS-] 



^{AB}S- 

-1{A,B}S+ 
{A,B}S^'+^-[A,B] 



(Bl) 



are fulfilled for arbitrary (linear) operators A, B that 
commute with the spin operators. 

Using Eq. (A7l the following relations are easily de- 



rived. For the 1-loop calculation, the commutator 



c\,ci 



4'C2 :] = : c\,C2 ■■ (5i,2'- : 4'^! ■ ^i'.2 (B2) 
+(5i,,2<5i,2'Wl')-"(l)) 



is needed. Due to the spin operator algebra (Bl ) also the 
anticommutator has to be calculated: 



{: 4'Ci : 4-^2 :} = 2 : 4'Ci4'C2 : 
-|-Ji,2'(l-2n(l)) :4,C2 : 

+(52,i'(l-2n(l')) :4'Ci : 

+<5i^2'5i,2' (n(l') + n{l) - 2n(l')n(l)) . 



and the anticommutator 

{: 4'Ci :, : 4'C24'C3 :} = (B5) 

<5l',25l,2'Wl')+^l) 
-'^l'.25l,3'Wl')+^(l) 

-Sv,3Si,2'in{l')+n{l) 



2n(l')n(l)) 
2n{l')n{l)) 
2n(l')n(l)) : c^,C2 
2n(l')n(l)) 



4'C3 



4'C2 



for the further calculation. 



Using the above relations the task of deriving the flow 
equations is reduced to simple but lengthy bookkeeping. 
The resulting 2-loop equations are given in the following. 
In the diagonal part of the Hamiltonian only the splitting 
of the dot levels due to the magnetic field is shifted 



dh 
dB 



(B6) 

x{e,~e, + h){J^„r + 0{J^) . 



In the case of zero (initial) magnetic field the relation 
Jpq = "^i^ is fulfilled leading to dh/dB — and therefore 
no additional magnetic field is generated. In the inter- 
action part we have to keep track of different scattering 
processes that lead to different flows of the running cou- 
plings though we started with isotropic initial conditions. 
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For the scattering of spin up electrons we find 
djT 1 

r 

- + "'/(^) ~ 2n/(r)n/(s)) J;L, ((e^ - + 2(er - e, + - Kl^,rq) 

-(ep - - 2(e, - e, + h)){Kl^^^ - kI^J) + 0{j') 
and for spin down scattering 

^ = -(ep - e,)V^ + - ^(2(e. + h) - (e^ + e,)) J.-^, ^,^^,(1 - 2n/(r)) (B8) 
+ Y.^rif{r) + nf{s) - 2nf{r)nf{s))J^, {{cp - + 2(e, - + - i^.^pJ 

-(ep - - 2{er - e, + - Z^^^,,,)) + 0{j') . 

The spin flip scattering is given by 

= -(ep - + + J 5^(1 - 2n/(r)) ((2e, - (e^ + e,) + /i) J^VJr + (2er - (e^ + e,) - /i) Jp"", J,V) (B9) 

r 

+ 1 Yi^fir) + nf{s) - 2nf{r)nf{s)){ep - + 2(e^ - e^) + h) 

r,s 

(("^pq,rs -^rq^ps)^ sr ^-^pq^rs ^ps,rq)'^sr} 

- J E^'^i^^) + "/(^) - 2n/(r)n/(s))(ep - + 2(e, - e«) - h)K^^^,,Ji + 0{J^) . 



""^pq 

~dB 



I 

The flow of the newly generated interactions is given by function and the response function. The flow can be 

described using one set of running couplings for both x- 



pq,TS 

dB 



-{cp -eg + er-es + hf^Klq ,^ 



(BIO) 



+ 4(ep - eg - er + e, + /l) JpVrs + ^("^^) 



for the spin up plus spin flip scattering and 



pq.rs 

dB 



-{Cp - Cq + Cr - eg + h)'^K^g . 



pq,rs 



(Bll) 



1 



+ e, + /i)J-44, + 0(j3) 



^q --r 1 1 ■"/•^pq'- 

for spin down plus spin flip. For double spin flip we flnd 



dK^ 

pq,TS 

dB 



'(Cp Eg + ^s) -f'^pq.rs 



-{ep -eq-er + es + 2h)J^qJ^^ + 0{J^) . 



APPENDIX C: TRANSFORMATION OF S"/'" 



and y-direction since the choice of the basis in the xy- 
plane is arbitrary. Also the correlation and the response 
function are given by single functions for both directions. 

We use the following ansatz for a;-direction 

S%B) = h^y{B)S^ + iY,f^lq{B) : fl^fq, : ^^Cl) 

+iY,niq{B):fljq^:Sy 

+ Y.Hlq{B){:fl,fq,: + :fyp, :)S^ 



(B12) and y-direction 



sy{B) = h-y{B)sy -iY^t^Us) f^fq, (C2) 

-iJ2f,i^{B):fyq^:S^ 

P.I 

-iY.t^lq{B){:fyq,:-:fyp, ■.)S^ . 



For completeness we show the transformation of the 

spin operators perpendicular to the magnetic field and 

give the result for the corresponding spin-spin correlation The flow equation for the decay of the spin operator is 
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given by 

^ = ^Y.Mp) + nf{q)-2nf{p)nf{q)) (C3) 

+ 1 + ~ 2n/(p)n/(g)) 

p,q 

The flow of the newly generated operators is given by 



pq 
dB 



2 'i^p ^'q)Jpq 



(C4) 



prl^qr 



■jE(l-2n/(r))(ep-e. + /i)Jp^, 

r 

4 11(1 - 2n/(r))(e, - + /i) J^i/i; 



-L ..z 
pr 



for spin down. For the spin-flip component we flnd 



pq 



hxy 



dB 



{ep - eg + h)J^g 



(C6) 



r 



T 



for the spin up component and 



pq 
dB 



"(Cp ^q)Jpq 



(C5) 



-L .,z 
rp 



- J - 2n/(r))(er - + /i) J^V 

r 

+ - ^^(1 — 2n/(r))(er — + h)J:^/j,^q The correlation function is given by the lengthy formula 



p 



p 

P 



4 

p 
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For the imaginary part of the response function we find 

7r(l + sgn(/i)) 



p 





8 


7r(l 


- sgn{h)) 


7r(l 


8 

- sgn(h)) 



p 

+ T 5I(^^P:Ep+'^)^("/('^p)(l - n/(ep + uj)) - nf{ep + w)(l - n/(ep))) 



4 



p 



Note that the equations above are identical to the transformation of the 5^ operator at zero magnetic field. 
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